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Abstract 



Deryagin, Grigoriev, and Rubakov (DGR) have shown that in finite-den- 
sity QCD at infinite the Fermi surface is unstable with respect to the 
formation of chiral waves with wavenumber twice the Fermi momentum, while 
the BCS instability is suppressed. We show here that at large, but finite A'c, 
the DGR instability only occurs in a finite window of chemical potentials 
from above Aqcd to /icrit ~ exp(7ln^A'^c + 0(lniVc Inln A^c))Aqcd5 where 
7 « 0.02173. Our analysis shows that, at least in the perturbative regime, 
the instability occurs only at extremely large A'c, ^ 1000 A'f, where A'f is 
the number of fiavors. We conclude that the DGR instability is not likely to 
occur in QCD with three colors, where the ground state at finite density is 
expected to be a color superconductor. We speculate on the possible structure 
of the ground state of finite-density QCD with very large A^c- 
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I. INTRODUCTION 



In contrast with finite-temperature QCD, QCD at high baryonic densities remains re- 
markably poorly understood. One of the main reasons is the lack of lattice simulations 
due to the complex fermion determinant in finite-density QCD. Meanwhile, the physics in 
the core of the neutron stars, and possibly of heavy-ion collisions, depends crucially on the 
structure and properties of the ground state of QCD at finite densities. 

It was suggested that at sufficiently high densities, the ground state of QCD is a color 
superconductor Such state arises from the instability of the Fermi surface under the 

formation of Cooper pairs of quarks. The superconducting phase of quark matter is the 
subject of many recent studies and we will not discuss its properties in this paper. We 
will only note that a reliable treatment is currently available only in the perturbative regime 
of asymptotically high densities 0]; in the physically most interesting regime of moderate 
densities, QCD is strongly coupled and one has to resort to various toy models, e.g. those 
with four- fermion interactions 0. 

To shed light on possible new phases that may occur in the non-perturbative regime of 
moderate baryonic densities, one might hope to be able to make use of alternative limits, 
such as the large Nc limit, where one takes the number of colors Nc to infinity, keeping g'^Nc 
fixed {g is the gauge coupling) [^]. This limit has proved to be a convenient framework for 
understanding many properties of QCD (for example, Zweig rule) , although QCD at infinite 
Nc is still not analytically treatable. In the context of finite-density QCD, the first work 
that discussed the implications of the large Nc limit was done by Deryagin, Grigoriev, and 
Rubakov (DGR) DGR noticed that color superconductivity is suppressed at large Nc 
due to the fact that the Cooper pair is not a color singlet (the diagram responsible for color 
superconductivity is non-planar) .0 Working in the perturbative regime g'^Nc <^ 1, DGR 
noticed another instability of the Fermi surface, this time with respect to the formation of 
chiral waves with wavenumber 2pp, where pp is the Fermi momentum. As shown by DGR, 
this instability is not suppressed in the limit iVc — >■ cx). 

The purpose of this paper is to see what happens to the DGR instability at large but 
finite Nc- Our motivation is to see whether the limit A^c oo is relevant for the physics of 
high-density QCD at A^c = 3. In this paper, we find that at any fixed value of the chemical 
potential fi, in order for the DGR instability to occur we require the number of colors Ac 
to be larger than some minimum value Ac(/i), which grows with /i. What is surprising is 
that even for moderate values of fi, the minimum value Ac(/i) is very large (of order of a few 
thousands for a modest chemical potential fi = 3Aqcd)- Therefore one should not expect 
the large Ac limit to be of direct relevance for physics with A'c = 3 at finite densities. 

The paper is organized as follows. Section |I| reviews the results of DGR. A convenient 
technical approach to DGR instability which is based on renormalization group is developed 
in Sec. ^TT| and applied to the case of finite A'c in Sec. Section |V| contains concluding 
remarks. 



-•^At arbitrary using the technique of Ref. |4|,the asymptotic behavior of the BCS gap can be found to 
be A ~ ^expl — y ^^^^ )• This tends to as A'c — > oo, provided one keeps g^Nc fixed. 



1 



II. REVIEW OF DGR RESULTS 



Let us review the key results of Ref. 0. Throughout our paper, we assume all quarks 
are massless, and make no distinction between Fermi momentum and Fermi energy: pp = ix. 
In the Nc ^ oo limit, the DGR result states that the Fermi surface is unstable under the 
development of chiral waves with wavenumber 2/i, 

(V^(a;)V(2/)) = e^P(^+y) 1^ e"^^^^"^V(g) (1) 

where P is a vector with modulus |P| = /i whose direction is fixed arbitrarily. Since ipip is 
a color singlet, it survives the limit N^. ^ oo. 




FIG. 1. The particle- hole pair 



The condensate (|I|) can be interpreted as the formation of particle-hole pairs with total 
momentum 2P (Fig. |I|). In such a pair, both the particle and the hole are near the Fermi 
surface, and the momenta of the particle and the hole are both near P. In this sense the 
condensate (0) is different from the usual chiral condensate (V'V') = const, which corresponds 
to the pairing of a particle and a antiparticle moving in opposite directions. Moving in the 
same directions, the scattering between the particle and the hole is nearly in the foward 
direction, and since the amplitude of forward scattering is singular, one could expect the 
formation of the pair to be energetically favorable. In fact, this is the reason why the total 
momentum 2/i is special. 

The function /(g) has the physical meaning of the wave function of the pair in the center- 
of-mass frame, so P + q is the momentum of the particle and P — q is that of the hole. 
DGR found that the wave function is localized in an exponentially small region of momenta 
q < A_L where 

A, ^e-"^\ = (2) 

Recall that h is kept constant in the limit ^ oo. The binding energy of the pair is found 
to be at an even smaller scale, 

^bind ~ /ie--/^ (3) 

Both scales A_|_ and -Ebind are parametrically larger than the non-perturbative scale Aqcd ~ 
^g-6/ii/i _ Pqj, ixiore details, see Ref. 0. 
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It may seem surprising that the DGR instabihty occurs in the perturbative regime. 
Indeed, the analogs of (|lD, in non-relativistic fermion systems, are the charge-density wave 
(CDW) and the spin-density wave (SDW). Since it is known that in three dimensions CDW 
and SDW do not develop at small four-fermion interaction, one could ask how such instability 
could occur at small g'^N^. The key observation is that in our case the effective four-fermion 
interaction is singular due to the behavior of the gluon propagator at small q. In 

general, this singularity is cut off by screening, but because the diagrams responsible for 
screening involve fermion loops, the screening effects at large are of order g"^^^ ~ 0{1/ N^^ 
and therefore suppressed. This singular nature of the interaction explains why the DGR 
instability can occur perturbatively at large N^. 

The argument presented above also implies that at each value of the coupling h, there 
must be a lower limit on N^., below which the interaction is not singular enough due to 
the screening, and the DGR instability disappears. This limit grows as one decreases h, 
or, equivalently, as the chemical potential increases. Finding this lower bound on as a 
function of /i is the purpose of this paper. 

III. RENORMALIZATION GROUP APPROACH TO DGR INSTABILITY 

Before tackling our main problem, let us formulate an efficient RG technique that repro- 
duces the results of DGR in the Nc ^ oo limit. While in the limit Nc —>■ oo this technique 
does not give us anything new over what has been already found by DGR, it has the advan- 
tage that it can be applied to the case of finite Nc, where the effects of screening make the 
generalization of the original method of Ref. ||^ very difficult, if at all possible. We will not 
try to rigorously justify the RG in this paper. 

Let us stay in the Fermi liquid phase, where quarks are deconfined, and consider the 
scattering between a particle and a hole with momenta P+q and P— q. The total momentum 
of the pair is 2/x. A singularity of this scattering amplitude in the upper half of the complex 
energy plane would signify an exponentially growing mode, i.e. an instability [0. In terms 
of the diagrams, the most important contribution to the scattering amplitude comes from 
the ladder graphs (Fig. |^). Adding a rung to the ladder brings two more logarithms: one 
comes from the coUinear divergence, i.e. the singular gluon propagator, and the other from 
the fact that the two new fermion propagators are near the mass shell. We will design the 
RG to resum these double logs.| 




FIG. 2. Ladder approximation 



similar but not identical RG procedure has been developed to resum the double logs in the BCS channel 

i- 
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Our first step is to derive a 1+1 dimensional effective theory capable of describing the 
DGR instability. On the most naive level, such description exists due to the fact that 
the modes of interest move in directions close to ±P. Technically, the (1+1)D effective 
theory arises from integrating, in each Feynman graph, over the momentum components 
perpendicular to P [Q. 

Let us consider a ladder graph and ask what happens if one adds one more rung. The 
diagram now contains an extra loop integral, 

j^,G{P + q)G{-P + q)D{q) (4) 

where G and D are fermion and gluon propagators respectively, and the Dirac structure 
of the fermion propagators is ignored for the purposes of the discussion presented below. 
Consider first the fermion line with momentum P + q. The component of the fermion 
momentum parallel to P will be denoted as /i + gy, and those perpendicular to P will be 
denoted as q±. Note that is a two-dimensional vector. We will assume that for all fermion 
lines in the Feynman diagram ~ A, where A is an arbitrary momentum scale much less 
than /i. In other words, we will be interested only in the modes located inside two small 
"patches" on the Fermi sphere, each having the size of order A in directions perpendicular 
to P (eventually, A will be identified with Aj^ in Eq. (^). When q\\ is also small compared 
to /i, the fermion propagator has the form 

G{q) r;^^ ^ ~ ^ . (5) 

ZgO+|P + q|-/i + + 

If q\\ ^ Q'i//^ ~ '^^//^j the g_L dependence drops out and the propagator is simply {iq^ + 
Therefore, in the regime q\\ ^ A^//i, the fermion propagator does not depend on 
the perpendicular (with respect to P) momenta. In this regime, in Eq. (^ only the gluon 
propagator D{q) depends on g^. Hence, the integration over q± has the form 

cPq^ 1 



(2vr)2 ql + ql + qV 

If go and q\\ are not only small compared to /x, but also much smaller than A, then the 
integral over g^ in Eq. (|^) is a logarithmic one / <Pq±/q\- The integral is cut off in the IR by 
q\\ and in the UV by A and yields |^ In ^. Effectively, this integration replaces the internal 

gluon line by a four- fermion vertex |^ In — , where gy is determined by the momentum of the 

fermions coming in and out of the vertex (Fig. Recall that the simplification takes place 
only in the region A ^ g|[ S> A^//i, as only in this region the integration over g_|_ decouples 
from that over gy. At the end of this section we argue why the restriction of gy to the region 
A^//i < gy < A is well justified. 
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We have taken the integration over the perpendicular components of one particular gluon 
momentum, but nothing prevents us from integrating over the perpendicular components of 
all the gluon momenta. By doing this integration, we resum one set of logarithms (the one 
related to the collinear divergence) in a series of double logs. Now, as the only remaining 
integrals are over go and q\\, all Feynman diagrams are identical to those of some 1+1 
dimensional model with a four-fermion interaction. Our task is to find out the precise form 
of the Lagrangian of this model. 

First, we note that the kinetic term for the fermions in the effective theory can be 
obtained from the original Lagrangian by omitting spatial derivatives in directions other 
than z, 



^kin 



It is more convenient, however, to recast the Lagrangian (|^) into the form of a (1+1)D 
theory of a doublet of Dirac fermions (which are two-component in (1 + 1)D) at zero chemical 
potential. This is indeed possible, since spinless fermions at finite chemical potential can 
be rewritten as one Dirac fermion at zero chemical potential (the modes near two points of 
the "Fermi surface" serve as its two components 0). It is not surprising that in our case 
the spin-| fermions can be rewritten as a doublet of (1+1)D Dirac fermions. Let us do it 
explicitly when P is directed along the z-axis, P = (0,0,/i). Denote the four components 
of the Dirac spinor ip (in chiral basis) as ip^ = (V'li, '?/'L2, ^ri, '?/'R2)- The antiparticles have 
energy of order 2/x and decouple from the low-energy effective theory that is being derived. 
This allows us to consider only the components of ijj corresponding to particles, which are 
'?/'l2 and ipYii when the particle's momentum is near P, and t/'li and '?/'r2 when it is near 
—P. Although these fields are slowly varying in time, they still vary rapidly in space. To 
compensate for this spatial variation, we introduce new fields, 

which are soft in both space and time. We can now translate from the (3+l)D language of 
ip to the (1+1)D language of ip and x- The kinetic part of the Lagrangian (|^) becomes 

^kin = iipl^doip + i^p'y^ds'ip + /iV'7°^ ^ i'^lm^f,^ + iXIm^f^X- (9) 

What is the interaction term in the effective theory? A look at the Feynman diagram in 
Fig. tells us that such interaction is of the current- current type. The current operator can 
also be translated into the (1 + 1)D counterparts, 

^720^^ + X72dX (10) 

where are two (1+1)D Dirac matrices, 720 = cr^, 720 = —ia"^. Below we will write 
these matrices simply as 7'^ in all expressions belonging to the (1+1)D effective theory. 
Noting that each vertex in Fig. |^ corresponds to a factor of |^ In ^, where gy is the parallel 

momentum transfer, we find that the Lagrangian of the (1+1)D effective theory is similar 
to that of the non-Abelian Thirring model 
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where we have combined the two fields Lp and x iii^o a doublet ^. The only difference 
between (|TT1) and the non-Abelian Thirring model is the dependence of the four-fermion 



coupling on the scale of the parallel momentum exchange gy . The theory ([TT|) describes the 
interaction between fermions with perpendicular momenta of order A and parallel momenta 
between j \x and A. 

To understand the properties of the model (|ll|); let us recall what is known about the 
conventional Thirring model, where the interaction term is — A(^7'^^^)^. The Thirring 
model is asymptotically free. The only diagram contributing to the /3 function at large 
is the "zero-sound" diagram, Fig. ^. 




FIG. 4. The zero-sound diagram 
The running of the coupling A is governed by the RG equation. 



ds TT 

where s is the RG parameter, and A(s) is the coupling at the energy scale Ae~^. The 
coupling A hits a Landau pole at p ~ /S.e'"^^^^" . The physics in the IR is characterized by 
the formation of the chiral condensate (\1'\E') which gives mass to the fermions. Using Eq. 
(I), one can see that (\E'\E') = cos 2nz{^p^) —i sin 2/iz(?/'7°7^?/'), so a constant (\E'\E') translates 
into space-dependent condensates (ipip) and (ip'y^'y^ip). 

These basic properties hold for the model ([TTl) as well, but the estimation for the scale 
of the Landau pole is different. The latter can be found using RG. Now the RG equation 
needs to be written for a coupling which is a function of the parallel momentum transfer gy. 
At s = 0, 

A 

A(g||) = ^ln-. (12) 
Air g|j 

The RG equation is found from the diagram drawn in Fig. ^. The internal fermion lines 
have the momentum of order Ae"'', which is much larger than the momentum of the external 
lines, therefore the momentum transfer at each vertex is Ae~*. The RG equation, therefore, 
is 

|-A(.,g||) = ^A2(s,Ae-^). 
as vr 

It is convenient to use the logarithmic parameter u, defined by gy = Ae~", and rewrite the 
RG equation as 

^Xis,u) = ^X'is,s). (13) 

as TT 
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The initial condition (|T2]) becomes 

A(0,u) = f-u. (14) 

4:71 

One should note that at the moment s of the RG evolution, all fermion modes with energy 
larger than Ae~'^ have been integrated out, therefore the function A(s, u) is defined only for 
u > s. The solution to Eq. (13) with the initial condition Eq. (|14|) is 



Hs,u) = ^J{s) + ^{u-s) (15) 
where /(s) satisfies the equation 

§-J{s) = e + fis) (16) 

and h"^ = g^Nc/^n'^. Solving Eq. ([l^) one finds /(s) = /i tan /is, which hits a Landau pole 
at s = sl = 7r/2/i. The corresponding scale is Ej^ = Ae"'^^'^^. Recall now that RG evolution 
occurs for A^//i ^ q\\ <^ A. From this condition one finds that the Landau pole can only 
be achieved if A^/^u ^ i?L or A ^ ne~^^'^^. Under this constraint, the maximal value of i?L 
is achieved when A ~ yue"'^/^^, at which Ei^ = A^//i ~ iie~'^f^. Thus, the estimation for 
the Landau pole scale Ei^ and for A coincide with the result found by DGR for the binding 
energy of the particle-hole pair and the size of the pair wave function, Eqs. (0,^. 

Now, it is easy to demonstrate why we were justified to consider only the region A^//x <S 
q\\ <^ A in the argument presented above. On one hand, when q\\ drops below the scale 
A^/yU, we cannot neglect the dependence of the fermion propagator on g^, which now acts 
as a cut off for the RG flow. Hence, for ^ A^/ yU, there is no RG flow in the effective 
(1+1)D theory and the Landau pole is never reached. On the other hand, when gy becomes 
comparable with (i.e. A), we cannot neglect q\\ dependence in the gluon propagator. One 
can estimate the effect of such dependence by noticing that the four-fermion coupling in the 
effective (1-|-1)D theory (|TT|) now reads 

A(gii) = 

and the RG equation (0) becomes 

One can see that for q\\ ^ A, the RG flow in the effective (H-l)D is completely negligible. 
Therefore, to flnd the DGR instability we can restrict the values of q\\ to lie between A^//i 
and A. 

Having reproduced the DGR results by our RG procedure, let us turn to the case of 
large, but flnite N^. 
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IV. DGR INSTABILITY AT FINITE A^c 



The RG technique described above can be very easily extended to the case of large but 
finite A^^c- The effect of finite N^. is to cut off the IR singularity of the gluon propagator at 
small momentum exchange by Thomas-Fermi screening and Landau damping. The electric 
propagator becomes (g^ + m^)~^, and the magnetic propagator becomes (g^ + im'^\qQ\/q)~^ 
10| , where m is the Thomas- Fermi screening scale of order g^. If the screening mass m is 
smaller than the scale of the Landau pole found in Sec. i.e. fie~'^^'^, then our previous 
calculations are not affected. However, if m > fxe^'^^'^, we need to modify the RG to take 
into account the screening. 

The screening affects the integration over perpendicular components of the gluon propa- 
gators: before these integrals were cut off by the parallel exchanged momentum gy, now it is 
cut off by the largest scale among gy and m in the case of electric gluons, and among q\\ and 
m^/^gy^^ in the case of magnetic gluons. The effective (1-|-1)D theory is now a Thirring-like 
model with different scale-dependent couplings for the electric and magnetic interactions, 

Leff = t^Yd,^ - Ao(g||)(^7°Y*)' + Y*)'- ^^^^ 



where 



Ao(g||) = |-ln ^ 



Air max(g||,m) 

Ai(g||) = |— In — Y73~- 

47r max(g||, m^/^gii ) 

The RG equations for A+ = (Aq + Ai)/2 and A_ = (Aq — Ai)/2 decouple: 

^X^(s,u) = ^\l{s,s) (19) 

OS TT 

where again u = In^. Therefore, only A+ changes during the RG evolution. The initial 
condition for A+ can be read from Eq. ([TsD , 

U it U < Sm 

M0.-)= p,, 1 . (20) 

where = In — . The solution to Eq. ( [T9| ) with the initial condition ( pOD can be written in 
the form of Eq. (|T5D , where /(s) now satisfies the equation 

Q ( p + h'^ if s < Sm 

d~s^^'^^[f + ^ its>sm ■ 
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The solution to this equation is 



h tanks if s < s 



m 



^ tan -^(s + c) if s > 
where c can be found by matching the solution at s = Sm- 



^ arctan(v^ tan /iSm) — 



The Landau pole occurs at 



V^TT -\/6 ,1 

sl = -^TT c = — arctan(^ cot hsm) + Sm- 

2n h v6 

Recall that for the instability to really occur, the scale of the Landau pole should be larger 
than the scale A^/ yU, one finds a condition on m, 



m = Ae """^ < fie ''^ = exp 



— — arctan —= cot hSm — 2sr 
h \^/6 



One can maximize the right hand side (RHS) of this equation to find the maximum value 
of m where the Landau pole still can be achieved. One finds that for the Landau pole to be 
reached, m should be smaller than mmax = fie~'^^'^, where 

c = V 6 arctan 2 ^ arctan W g ~ 2.5051. 

This restriction on m leads to a condition on A'^c and fi for the DGR instability to occur. 
Recall that the Thomas-Fermi mass is 



27^ 



(which is of order Nc'^^"^), we see that at a fixed coupling g'^Nc (or, equivalently, /i), there 
exists a lower bound on Nc where condition m < /le"'^/'^ is satisfied. The lower bound can 
be easily found to be 

iVe > 2Nih^e^^/''. (21) 

Since our arguments rely on the comparison of scales, Eq. ( pT| ) contains an extra unknown 
coefficient of order 1 on the RHS. As the chemical potential fi increases, the effective coupling 
h decreases; using the one-loop beta function 

= , ^ „ (22) 

llln-r^ 

Aqcd 

and according to Eq. (^) the lower bound on A"c increases. In reality, the numerical constant 
2c in the exponent on the RHS of Eq. (|2T|) is relatively large (~ 5), so the lower bound is 
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already large at moderate values of /i. For example, if one uses the value of h corresponding 
to /i = 3Aqcd, the RHS of Eq. is of order lOOOA^f! Barring the possibility of a very 
small numerical constant on the RHS of Eq. (0), which seems unlikely, this lower bound is 
always much larger than 3. 
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ln(n/AQ^j^) = y ln^(1 1/12 N_,/N, ln(n/A^ 




2000 4000 6000 8000 10000 12000 14000 16000 18000 



FIG. 5. Region of DGR instability in the {Nc,fi) plane from Eq. (^). CS and DGR stand for regions 
with predominant color superconductivity and DGR instability, respectively. 



From Eqs. ( pl| , p2D , one can construct the phase diagram of QCD in the {Nc, jj,) plane. 
The result is shown in Fig. |^. In the shaded region, satisfies the inequality (pi]), which 
means that DGR instability occurs. We restrict this region by the line /z = 3Aqcd, for below 
this line QCD is certainly strongly-coupled and not much can be said from our calculation. 
Above the curved line, inequality (^TJ) is not satisfied, and the Fermi surface is stable in the 
DGR channel. However, the BCS instability is still there (though suppressed by large Nc), 
thus implying that the ground state of QCD is a color superconductor in that region. 

At any given (large) N^, the DGR instability occurs only in a finite window of the values 
of the chemical potential. The maximal value of fi where DGR instability still occurs, /^crit, 
can be found by solving (|21| ) with respect to /i. Asymptotically, 



/icrit ~ exp(7 In^ iVc + 0(ln In In A^c))Aqcd 



iVc^'^'^^AQCD 



where 
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22c2 



0.02173. 
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The smallness of the numerical constant 7 and the logarithmic dependence of /icrit on A'^c 
are the reasons why it requires a numerically large N^. for ^crit to be as small as 3Aqcd- 
However, asymptotically /icnt grows faster than any power of A*",,. 

V. CONCLUSION 

In this paper we have seen that in finite-density QCD the Fermi surface is unstable under 
the DGR instability in a finite range of chemical potential. We have also found that the 
number of colors needs to be numerically very large for the DGR instability to occur 
in perturbation theory. This indicates that at low Nc (like Nc = 3), the DGR instability 
might not have a chance to realize itself at any value of the chemical potential and the only 
instability of the Fermi surface is the BCS one, which leads to color superconductivity. 

Returning to the case of very large N^, the next logical step is to ask what is the ground 
state once the Fermi liquid is unstable under the DGR particle-hole pairing. This seems to 
be a purely academic exercise due to the large A^'c required, but it might still be interesting 
because of the possibility, at least in principle, of a new phase, distinct from the Fermi liquid 
and BCS superconducting phases in 3D fermionic systems. In the original paper 0, DGR 
constructed a "standing chiral wave state", in which {^tp) varies periodically in space with 
wavenumber 2/i. This state is periodic only along one spatial direction and does not break 
translational symmetry along the other two directions. Since translational symmetry cannot 
be broken in only one direction, such state cannot be the ground state of QCD. 

One notices that a chiral wave with a particular wavevector utilizes only fermion modes 
in a small region with size (Eq. (^) near two opposite points on the Fermi sphere. It is 
clear how to make a state with energy smaller than the original DGR standing wave state. 
Indeed, one can pair up particles and holes in different pairs of opposite patches on the 
Fermi sphere. Since the size of each patch is exponentially small compared to the total area 
of the Fermi surface, one can have a large number of patches that do not overlap with each 
other. From the size of the patches one deduces that one can place a maximum of e~'^^'^ 
patches on the sphere. The condensate has the form of a linear combination of e*''*'^, where 
all kj have modulus equal to 2/i but point in different directions. It is easy to estimate the 
energy gain from forming such a state. Indeed, the pairing affects fermions in a thin shell 
near the Fermi surface; the thickness of the shell is the scale at which we have found the 
Landau pole, i.e. /ie"'^/'^. Therefore, the fraction of fermions affected is e"'^/'*, and each pair 
lowers the energy by fie^^^'^. Therefore, the gain in energy density is 

/i^e-2-/^ (23) 

For comparison, the DGR standing wave state has the energy gain fj.^e~^'^^^. The factor of 
Q-->r/h flifference is explained by the fact that DGR state involves only two patches on the 
Fermi surface with a relative area of e'"^^^. 

Alternatively, it might be energetically more favorable for the patches on the Fermi 
sphere to be overlapping. In this case, a given particle (or hole) near the Fermi sphere 
participates in many pairings simultaneously. It could be expected that the binding energy 
of each individual pair is lower than the value it would have in the non-overlapping case, but 
nothing can be said about the total energy of the system. Indeed, our preliminary estimation 
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shows that the energy gain is still parametrically given by Eq. (]23|). Further investigation is 
required to find the true ground state of QCD at very large A^c- 

Finally, let us note an interesting possibility that the ground state of finite-density QCD 
at very large might be similar to the "tomographic Luttinger liquid" in 2D, advocated 
by Anderson as the normal state of high-Tc cuprates ||TT|. Such similarity could stem from 
the singular interaction between fermions moving in the same directions, which is also char- 
acteristic of tomographic Luttinger liquids. As in the case of the latter, one could expect 
the chiral symmetry to be unbroken, but the chiral response to be singular at wavenumber 
2/i. 
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